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A recent theory by Feigel [Phys. Rev. Lett. 92, 020404 (2004)] predicts the finite 
transfer of momentum from the quantum vacuum to a fluid placed in strong per- 
pendicular electric and magnetic fields. The momentum transfer arises because of 
the optically anisotropic magnetoelectric response induced in the fluid by the fields. 
After summarising Feigel's original assumptions and derivation (corrected of trivial 
mistakes), we rederive the same result by a simpler route, validating Feigel's semi- 
classical approach. We then derive the stress exerted by the vacuum on the fluid 
which, if the Feigel hypothesis is correct, should induce a Poiseuille flow in a tube 
with maximum speed « 100/im/s (2000 times larger than Feigel's original predic- 
tion) . An experiment is suggested to test this prediction for an organometallic fluid 
in a tube passing through the bore of a high strength magnet. The predicted flow 
can be measured directly by tracking microscopy or indirectly by measuring the 
flow rate (w lml/min) corresponding to the Poiseuille flow. A second experiment 
is also proposed whereby a 'vacuum radiometer' is used to test a recent prediction 
that the net force on a magnetoelectric slab in the vacuum should be zero. 

Keywords: Feigel effect, quantum vacuum, magnetoelectric, strong magnetic 
fields, particle tracking velocimetry, vacuum radiometer 



1. Introduction 

It is well known that quantum vacuum fluctuations can transfer momentum to 
macroscopic matter. The momentum transfer usually results from a modification 
the spectrum of allowed vacuum modes by symmetric boundaries. For example, in 
the Casimir effect (Casimir 1948; Lamoreaux 1997, 2005; Milonni 1994) two con- 
ducting plates are placed parallel to each other. The plated reduce of the allowed 
vacuum modes between them, which causes the total pressure (due to the momen- 
tum transfer of vacuum modes) to be smaller between the plates than outside, 
resulting in a net attraction between the plates: the Casimir force (Milonni et al 
1988). 

However, the vacuum does not usually transfer its momentum to matter in the 
absence of boundaries: because its fluctuation spectrum is isotropic in free space the 
expectation value of the momentum density of the vacuum is zero (Milonni 1994). In 
this paper we revisit a theoretical argument proposed recently (Feigel 2004), which 
claims that it is possible to transfer momentum to an isolated region of dielectric 
liquid if it is placed in perpendicular crossed electric and magnetic fields. Critical 
to Feigel's argument is the assumption that vacuum fluctuations see the dielectric 
in crossed fields as an anisotropic magnetoelectric medium. This implies vacuum 
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modes will propagate faster (have greater momentum) in one direction than in the 
opposite direction, causing a net momentum transfer to the dielectric liquid. 

In going through Feigel's original argument we have corrected trivial errors in 
the original derivation. The corrected expressions agree with our own more direct 
derivation of Feigel's result. In the original paper Feigel proposes a prediction for 
the speed of a fluid as a result of the momentum transfer from the vacuum. The 
prediction was, however, predicated on an ideal fluid, and thus not realistically 
testable. We derive here the flow of a real, viscous fluid driven into motion by the 
stress on the fluid caused by the vacuum. An experimental test of this improved 
prediction is proposed. Indeed, this paper is written in the hope that our revised 
predictions will stimulate experiments to try and measure Feigel's original effect 
and related recently proposed predictions. Since the publication of Feigel's result, 
several theoretical works have emerged which have been very enthusiastic about 
Feigel's idea (van Tiggelen et al 2005, van Tiggelen et al 2006, Shen et al 2006, 
Birkcland & Brevik 2007, Obukhov & Hehl 2008). Some have, however, questioned 
the soundness of the original argument stating that no unbounded macroscopic 
Feigel effect should exist if proper regularisation is applied to the momentum inte- 
gral (van Tiggelen et al 2005, van Tiggelen et al 2006). Microscopically, however, 
recent work suggests that a properly regularised Feigel effect could exist (Kawka 
et al 2010). Most other alternative macroscopic theories consider the momentum 
from the vacuum for magnetoelectric fluids (Birkeland & Brevik 2007) and samples 
of magnetoelectric materials (van Tiggelen et al 2006) confined in parallel plate 
geometries, like in the Casimir effect. The latter work predicts an unmeasurably 
small linear momentum transfer (van Tiggelen et al 2006). Interestingly, follow- 
ing a semi-classical approach similar to Feigel, Obukhov & Hehl (2008) found that 
there is no net force on a magnetoelectric slab of finite thickness in the vacuum. 
Feigel (2009) has also recently considered the interesting possibility of constructing 
'quantum wheels' using magnetoelectric nanoparticles. No experiments have been 
carried out to test whether Obukhov & Hehl's null result, which implies Feigel's 
wheels should not function. More surprisingly, Feigel's original theory, which re- 
mains the one predicting the largest effect, remains untested. Experimental tests 
for the effects described are suggested at the end of the paper, but we leave the 
field open to imaginative experimentalists with access to the strong fields or mag- 
netoelectric materials required. 

The paper is organised as follows: in Section [2] the physical set-up for the model 
is presented together with quick summary of the derivation and underlying assump- 
tions. In Section [3] our quicker route to Feigel's quantum result is then presented; 
in Section 0] a new expression for vacuum stress is derived and a revised predic- 
tion of Feigel's theory is evaluated and discussed together with Obukhov & Hehl's 
null prediction; in Section [5] we make some suggestions for experimentally testing 
Feigel's and Obukhov & Hehl's predictions. Finally, in Section [5] conclusions are 
drawn from the preceding analysis and a few final comments are made. 



2. Feigel's semi-classical model 

We summarise here Feigel's semi-classical derivation and assumptions (highlighted 
in italics). Slight inaccuracies in the derivation of two key results have been cor- 
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rected, so these expressions differ from those in Feigel's original paper (a table 
comparing the original results with our corrected ones is presented in |Appendix C[ ). 

(a) Physical system and initial conditions 

Feigel considers the following situation: a region of a dielectric fluid far from 
the boundaries of its container is initially at rest it = 0). Subsequently strong 
electric and magnetic fields crossed at right angles to each other are applied to the 
region. As the fields reach their constant final values, E cxt and B ext for electric and 
magnetic fields respectively, the fluid is accelerated by Lorentz forces (-Fhorcntz oc 
<9 f (E cxt x B oxt )) to a final velocity v. 



(b) Brief summary of Feigel's derivation and fundamental assumptions 

1. Feigel assumes that the portion of dielectric fluid under consideration is to 
a good approximation ideal (inviscid), incompressible, homogeneous and not 
acted upon by external stresses or body forces] 

2. as the fields become steady, they have a maximum momentum, the oppo- 
site of which is transferred to the fluid by momentum conservation, since the 
combined system conserves its initial zero net momentum. Feigel derives the 
conservation law from the relativistically transformed Lagrangian of the mov- 
ing dielectric. This provides the 'classical' fluid momentum (ignoring terms of 
order (v/c) 2 ) 

pv = ^l (E xB) (2.1) 

47T/iC 

It should be noted that this conservation law holds only if the region under 
consideration is stress free, as mentioned in point 1. Otherwise momentum is 
not conserved. 

3. It has been shown that the optical response of a dielectric in crossed E and 
B fields is the same as that of a magnetoelectric material (Roth & Rikken 
2002). Feigel assumes that electromagnetic modes of the vacuum will also 
'see ' a magnetoelectric. The Lagrangian density of a magnetoelectric mate- 
rial is derived by Feigel by relativistic transformation of the Lagrangian of a 
magnetoelectric (see equation (|B 7[) of | Appendix B I in the small speed limit 



(ignoring, as above, terms of order (v/c) 2 ~). The Euler-Lagrange equation of 
this Lagrangian then provides a momentum conservation law: 

pv = J_ ( € J!—Le x B + -E x ( X T E) + - ( X B) x b) (2.2) 

47T \ fiC flC /J.C ) 

which allows to evaluate the momentum associated with the vacuum. We 
note, that using the inverted constitutive equations (IB 31) and (IB4[) . we can 
rewrite the above as: 

pv = -J— (D x B - E x H) (2.3) 
Aire 

we will consider the appropriateness of using this momentum density to eval- 
uate the fluid momentum in the Discussion. 
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4. The field modes and refractive indices for electromagnetic waves in a mag- 
netoelectric are given by (taking the optical axis along the e3 direction, see 
|Appoiidix A I , 

[E±ki,B ±k i] = i?±ki[ei,n ± kie 2 ], [E ±k 2,B ±k 2] = -E±k2[e2, -n±k2ei] 
where E±\ IL \ — -Eofee^* 2- "'- 1 , for each polarisation A = 1, 2, and 

n±k,l = ±"0 + Xxy, ™±k,2 = ±n - Xyx 

where \ xy and \yx are the magnetoelectric susceptibilities responsible for the 
optical anisotropy of the magnetoelectric (see | Appendix A[ ). Feigel substitutes 
the above modes and refractive indices into (|2.2D to evaluate the time-averaged 
(denoted by the overbar) momentum flux for a mode in the z-direction: 

pv k = 2A X - e -^ (2.4) 

C 47T 

where Ax = Xxy ~ Xyx- 

5. Feigel then replaces the electric field with its operators and evaluates the 
electromagnetic vacuum energy density expectation value 

<o|#|o> = ^ 

to obtain the vacuum momentum density per mode from (|2.4j) : 

g Qk = (0\pv k \0) = ^A X hk Ol (2.5) 

where fcO = lj/c. 

6. Summing over all modes the total momentum density in the z-direction is 
then: 

1 1 f°° 

where the last step involves the standard replacement ^2 ko — > g^r / al 3 ko (see 
next section). 

7. The integral in (|2.6p is divergent. Feigel makes the crucial assumption that 
vacuum modes with frequency greater than the dielectric's "cutoff frequency", 
lo c ( the frequency above which the dielectric 's molecular polaris ability van- 
ishes) do not interact with it. (Implicit in the derivation is also the assump- 
tion that absorption and dispersion are not significant, i.e. that for uj < u c , 
the permittivity and magnetoelectric susceptibility of the dielectric as seen 
by vacuum modes does not change appreciably with frequency). This allows 
to evaluate a finite value for the momentum density 



1 Tiac . . 

g a = pv v&c = —A x — . (2.7) 
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From which, dividing by the fluid density p, Feigel obtained an estimate for 
the vacuum contribution to the fluid speed, v vac . We iterate that the above 
expression is different from that obtained by Feigel due to some trivial errors 
in the original derivation (see | Appendix B[ ) . 

8. The magnitude of the corresponding "classical" contribution of the dielectric 
fluid's speed is given by (12 . 1 1) 

1 e/j, — 1 

^class = ; ^ cxt -D cxt- (2-8) 

p 47T/IC 

the relative magnitude of quantum vacuum and classical contributions will be 
discussed in section |U 



3. Simpler derivation of the vacuum contribution to 
magnetoelectric momentum density 

Here we present our alternative route to Feigel's result, equivalent to steps f — 5, 
but very much quicker. From quantum electrodynamics the expectation value for 
the momentum density g of the vacuum is given by (Milonni f994): 

go^(0|g|0) = i^iftk, (3.1) 

kA 

where V is a sample volume of the medium under consideration and k is the wave 
vector of each vacuum mode in this medium with possible polarisation states A = 
1,2. Next we assume with Feigel that the vacuum experiences the same birefringence 
in a magnetoelectric medium as light does. The medium parallel to the optical axis 
in the e z direction has the following dispersion relation: 

k • e z = fc «kA, (3.2) 

where rikA (given below) are the refractive indices parallel to the optical axis and 
we define fco = lo/c. On the other hand, in directions perpendicular to the optical 
axis the medium is isotropic, so that contributions to (|3.ip vanish by symmetry. 
Substituting in p.2[) then reduces to: 

go = y-fi'J2'J2ko(n +k x + n- k x)e z , (3.3) 

ko A 

where the sum has been expanded in terms of the contributions by modes for each 
direction (±) of travel along e z . The expressions for the anisotropic indices n±k\, 
derived in | Appendix A| are: 

n ±k ,i = ±n + Xxy, n±k,2 = ±n - Xyx, (3.4) 

where we recall \xij\ "C |no| so, e.g., rik,i > and n_k.i < 0. Substituting these 
expressions into (|3. 31) and summing over all polarisations A gives (considering only 
the magnitude of go): 

.9o = yAxh^2 k o, (3-5) 
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where A% = Xxy — Xyx- Next, making the standard replacement ^2 k — > J d 3 ko, 
we find: 



9o = A A ^ ft / k o dk °- ( 3 - 6 ) 



2vr 2 



i) 



These integration limits assume modes of all wavelengths contribute to the mo- 
mentum density, as they would in free space, so go is divergent. However, vacuum 
electromagnetic modes with very small wavelengths are not expected to interact 
with the macroscopic electromagnetic properties of a material medium. The choice 
of a reasonable value for the cut-off will be discussed in the next section. Here, 
with Feigel, we simply assume it is reasonable to approximate (|3.6p using an upper 
cut-off on the wavenumber k c — lj c / c to finally obtain: 

Equation (13.71) represents the vacuum contribution to the momentum density of a 
magnetoelectric. In media without the special symmetry of magnetoelectrics the 
refractive indices obey n+kA = — w>— k\, so that = and there is no transfer 
of momentum, as expected. Our alternative derivation provides a concise route to 
Feigel's result. It is equivalent to his semi-classical approach (once trivial errors in 
the original derivation are corrected, see | Appendix C| . 

4. Vacuum stress and realistic predictions 

The original prediction of the Feigel effect was obtained assuming conservation of 
momentum to obtain the vacuum contribution to the speed of a dielectric fluid 
placed in crossed fields (an effective magnetoelectric) from (|3.7p (Feigel 2004). 
Feigel's fluid was ideal and far from boundaries, making the measurement of the 
vacuum speed as originally predicted a Utopian pursuit. We show here how more 
realistic predictions for the original experiments can simply be obtained and pro- 
pose a new experimental test. These improved predictions are based on the fact 
that the transfer of momentum from the vacuum to a magnetoelectric results in a 
stress. 

For the case of an effective magnetoelectric fluid, an expression for the vacuum 
stress can be derived applying kinetic theory to a gas of vacuum modes (virtual 
photons) of momentum ^?ik travelling in the fluid. Optical anisotropy implies that 
the net momentum 

APk = 2 hk o^2 ("+ kA + n -kA)e 3 (4.1) 

A 

is transferred by counterpropagating vacuum modes across a surface in the fluid of 
area A and normal to e z . The corresponding stress on the fluid is the time rate of 
change of this momentum transfer per unit area: Ilk = \ ■ Modes crossing A in 

the interval At are recruited from a slice of fluid of thickness Az — c\ — — ]At 

(recall n_kA < 0), where c is the speed of light in vacuo. Thus, the magnitude of 
the net stress in the z-direction due to a mode pair is given by: 



n k = --^fc c> (4.2) 

V 2 ^ n + kA^-kA 
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where V — AAz. Substituting the refractive indices p. 41) into equation (I4.2j) (ne- 
glecting terms of order ||xl| 2 ) and summing over mode magnitudes gives: 

IIo = 2 — ^AxftVfco. (4.3) 
n V ^ 

fc 

Comparing this expression with (|3.5I) we see that IIo = 2goc/no, as one would 
expect (Loudon et al 2005). We can then use (|3.7[) to write: 

n = 4^ 2 A X ft— 4- 

where we recall no = ^Jep is the fluid's isotropic refractive index and note that (|4.4[) 
has been re-written in terms of the cut-off wavelength A c in view of the evaluations 
here below. Equation (|4.4p is the magnitude of the stress (acting in the z-direction) 
exerted by the vacuum on an effectively magnetoelectric fluid. In section (a) below 
we apply this expression to calculate the speed of a dielectric fluid in a tube with 
a portion of its length placed in strong crossed fields (a realistic version of Feigel's 
original scenario). Note that an equivalent result for the vacuum stress could also 
have been obtained semi-classically by evaluating the contribution from counter- 
propagating modes to the electromagnetic stress tensor: H^ z ~ eEofa-icA — n%. x ). 

Indeed, it was by using such a semi-classical approach applied to a magneto- 
electric slab of finite thickness in a vacuum that Obukhov & Hehl (2008) recently 
predicted that the net stress on the slab should be 

n = 0. (4.5) 

That is, the vacuum exerts no net force on a magnetoelectric the slab in spite 
of its anisotropy! In (b) below we propose to test this prediction with an experi- 
ment where the angular drift speed of a 'vacuum radiometer' with paddles made of 
magnetoelectric materials is measured. 



(a) Dielectric fluid in a tube in crossed fields 

We consider Feigel's original situation of a dielectric fluid placed in perpen- 
dicular electric and magnetic fields. In our case, however, the fluid is realistically 
contained in a tube, a section of which is exposed to the field (see figure [T^l). When 
the fields are switched on, Lorentz and ponderomotive forces act on the fluid. The 
fields induce a magnetoelectric susceptibility in the fluid where the fields act, so, 
according to Feigel's theory, the vacuum exerts a stress in that region. The fluid in 
the tube thus obeys the following Navier-Stokes equation: 

dv dP 
PM— = (P- V)E+— x B + 77V 2 v + V-n , (4.6) 

where v is the flow speed, E and B are the imposed fields, P = (e — 1)E + xB/// + 
°(llxl| 2 ) is the electrical polarisation of the fluid, pu and 77 are the density and 
kinematic viscosity of the fluid, respectively, and IIo is the stress on the fluid due 
to the vacuum. The fluid is assumed incompressible so V • v = and we note that 
the advection term (v • V)v vanishes in the cylindrical geometry of a tube. When 
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Planck's angular constant 


h 


1.05 x 10" 34 Js 


speed of light 


c 


2.99 x 10 s ms" 1 


permittivity of free space 




8.85 x tCP^Fm -1 


external electric field 


-Eext 


10 5 Vm" 1 


external magnetic field 


Bcxt 


17T 


ME susceptibility diff. (632.8 nm, room T) 


Ax 


lCT 11 (Roth & Rikken 2002) 


kinematic viscosity (25° C) 


V 


4.5 x 10~ 3 Pas (YiXing-KaiRun 2011) 


density (25° C) 


PM 


1.38 gcm~ 3 (Sigma) 


refractive index (589.3 nm, 20°C) 


no 


1.58 (Sigma) 


cut-off wavelength (cis-Pt-DEBP) 


A c 


4 nm (Matassa et al 2010) 



Table 1. Values of the parameters used in the evaluation of the vacuum and classical 
contributions to the velocity of an effectively-magnetoelectric organometallic fluid. 



both fields and the flow are in the steady-state, and ignoring ponderomotive forces 
due to edge effects, equation (|4.6p reduces to nV 2 v + V • n = 0. We thus expect a 
standard Poiseuille flow solution, with a maximum speed, U max = max|v|, at the 
center of the pipe given by (Brody et al. 1996): 

,„„ - ,,7, 

where IIo is the magnitude of the vacuum stress, a and L are the tube diameter 
and length, respectively, and r\ is the dynamic viscosity of the fluid, as above. 

To estimate U max we consider, like Feigel, the same organometallic liquids and 
the same magnitudes of crossed E and B fields used in the experiments by Roth and 
Rikken (Roth & Rikken 2002). In particular we focus on methylcyclopentadienyl 
manganese tricarbonyl (MMT), whose relevant properties are shown in TableQ] The 
magnitude of the vacuum stress is given by (|4.4[) and is seen to depend sensitively 
on the cut-off wavelength A c . Feigel used A c ~ O.lnm, taking intermolecular dis- 
tances as a cut-off. However, at such a scales vacuum modes will not interact with a 
sufficiently large number of molecules to experience magnetoelectric anisotropy. For 
larger distances, on other hand, modes can interact with the collective effect of an 
assembly of molecules. The radial distribution function (RDF) for a fluid provides 
a good measure of the distance beyond which modes see a smooth electromagnetic 
landscape. The RDF for MMT has not been measured to the best of our knowl- 
edge but a recent X-ray diffraction study has measured it for the powder similar 
organometallic compound, cis-Pt-DEBP (Matassa et al 2010). In the study it is 
found that the RDF becomes flat (uniform density) for > 4nm, which we assume 
as the value of the cut-off A c . The difference in susceptibilities A% is estimated 
approximately from the value of the birefringence measured by Roth and Rikken 
experiment (Roth & Rikken 2002). Using Equations (I3.4|) . we see that this bire- 
fringence is An = Xxy + Xyx which Roth and Rikken measured as An « 10~ n for 
the applied fields shown in Table [1] With Feigel we approximate A% w An. Hence, 
using the parameters of Table [T] equation (14.4)) gives a vacuum stress of Ho = 0.03 
Pa, which using (I4.7j) . for a tube with a — 1 mm and L ~ 2 m in (|4.7p . implies a 
flow with maximum speed 

U max = lOO^m/s- (4.8) 
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This prediction for the flow speed is 4000 times larger than Feigel's original predic- 
tion. The corresponding "classical" contribution to the velocity is negligibly smaller, 
7.5 nm/s, as can be seen from (|2.8|) using the parameters in Table [TJ As a matter 
of interest we note that (12.81) follows from (|4.6p when ponderomotive, vacuum and 
viscous stress contributions are neglected, so that the classical momentum is con- 
served. Inclusion of these realistic contributions, however, shows that the classical 
contribution, which we recall is due to Lorentz forces caused by polarisation cur- 
rents, is a transient, so it is not only negligible in magnitude, but also at long times 
when fields are steady. It will not contribute to the steady flow predicted above. 



(b) A vacuum radiometer 

Another way to test the reality of vacuum momentum might be to use natu- 
rally magnetoelectric compounds. In the organometallic fluids just discussed it is 
the high fields which induce magnetoelectric anisotropy in the fluid microstruc- 
ture. However, this anisotropy, and the resulting birefringence, can also arise in 
solids whose structures support both spontaneous polarisation and magnetisation, 
breaking both time and space inversion symmetries (Figotin & Vitebsky 2001). 
By virtue of this birefringence a slab of such materials in the vacuum should ac- 
quire momentum from the vacuum. Now consider an arrangement analogous to a 
Crookes' radiometer (Crookes 1876, Woodruff 1968), but with the vacuum pressure 
driving rotation as opposed to temperature gradients: a mill consisting of two or 
four square panes made of thin magnetoelectric slabs joined together with rods 
hinged on a low friction axle (see Fig. QJi). If each pane has area A and distance I 
from the centre of rotation, a vacuum stress of magnitude IIo normal to each pane 
causes a torque To = HqAI. The equation of motion of the radiometer mill is then: 

Id) = —^juj + tq. (4-9) 

where 7 is the frictional damping constant and I = p s A8l 2 its the moment of inertia, 
where p s is the slab density and <5 its thickness. Integrating equation (|4.9I) . we find: 

u(t) = Woo (l-e-*/ tc ). (4.10) 

where Woo = T0/7 is the terminal angular speed and t c = 7/7 is the characteristic 
time for the approach to this speed. Assuming the prediction of Obukhov & Hehl 
(2008) also applies to slabs of finite extent (i.e. the contributions of any stresses at 
the edge of the slab are negligible), we then expect 

u(t)=0. (4.11) 

The radiometer should not turn. If in an experiment it did actually turn, measuring 
its angular drift (see below) would allow to estimate the vacuum stress from 

n = • (4.12) 

We would expect such a stress should to scale with A%, as in (|4.4I) . 
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5. Possible experimental tests 

If the Feigel effect has the magnitude we have calculated, it should be possible 
to test the prediction for a dielectric fluid's velocity using current experimental 
techniques. Very recently, a 17 Tesla magnet with temperature control in the range 
1.6 — 300K and conical (±10°) access to its bore has been built (Holmes et al 2010). 
The magnet was designed for small angle neutron or X-ray scattering experiments, 
but it could also used to test our prediction (|4.8I) . Organometallic fluids would be 
placed in a tube arrangement going in and out of the magnet; a short portion of the 
tube would reside inside the magnet bore and would be fitted with electrodes to 
generate the required 10 5 Vm _1 electric fields. A schematic of the set-up is shown in 
Fig. QJ,. As well as electrodes which generate a large electric field, the high strength 
magnet would have to have a bore large enough to accommodate a tube for the fluid 
to flow in. In addition, fields need to be uniform over a channel and kept gradient- 
free to a good approximation. The 17T magnet described has 0.1% uniformity over 
lcm and a bore diameter < 4 cm (A. Holmes, Private Communication), so it would 
be reasonable to use tubes with ~ mm diameters, as assumed in scenario (a). Our 



(a) 



\ N / 














(b) 




Figure 1. Schematic of the proposed experimental arrangements for testing the Feigel 
effect, (a) A tube is placed in the magnet bore where it is subjected to a strong magnetic 
field and a perpendicular electric field. The channel enters and exits the magnet via a 
conical aperture, and flow is measured down stream by the elution or tracking the drift 
of colloidal tracers. Reversing the electric field allows to reverse the flow and eliminate 
unwanted drifts, (b) A 'vacuum radiometer' with magnetoelectric panes should rotate is 
the vacuum exerts a pressure on them. According to a prediction by Obukhov & Hehl 
(2008) this should not happen. An LED laser on the top of the slab allows to track any 
displacement of the radiometer in a darkened container (to avoid it being set in motion 
by ambient light. 



prediction is that the vacuum should generate a Poiseuille flow with maximum 
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speed 100/im/s, as obtained from (|4.7p . The corresponding flow rate, $ is given by 

$ = w 1ml mm , (5.1) 

so that, supposing the fields are as shown in Fig. [T^,, experimenters should be 
able to collect 1ml of fluid a minute after opening a tap at B (once the fields are 
in steady state). Reversing one of the fields, the same amount should instead be 
accumulate at B (this should also allow to subtract out any systematic variations). 
MMT is transparent, so, if the estimate (|4.7I) is incorrect but a finite vacuum 
induced flow still exists, it might also be possible to place stripped down (and 
minimally magnetic!) microscope with a camera down stream to perform particle 
velocimetry on colloidal tracers. The small flows inside the tube could then be 
measured downstream far enough away from the magnet to avoid perturbations by 
large stray fields on the velocimetry apparatus. The mean, (X), and mean squared, 
(X 2 ), displacements of these tracers could measured as a function of time, T: 

(X) = U T (5.2) 
(X 2 ) = ADT + U'lT 2 (5.3) 

Integrating measurements over a long enough time and fitting (|5.2[) and (|5.3[) to 
the data should allow to establish whether there is a nonzero vacuum drift Uq in 
the fluid and what its magnitude may be. From (15.31) . any vacuum induced drift 
present can be distinguished beyond doubt from Brownian motion by tracking for 
T 3> 4D/Uq. Supposing tracking can be carried out for for T exp , then, if no drift 
is observed, we infer: < Uq < AD /T exp . Again, the direction of the flow should 
be reversed by reversing one of the fields to eliminate other unwanted systematic 
drifts. 

Turning now to the prediction of Obukhov & Hehl (2008), this implies that a 
vacuum radiometer, such as the four pane one shown in Fig. [TJd, should not turn. 
An experimental test would involve a radiometer with thin panes made out of a 
solid magnetoelectric material (see Fig. [Tb). A good candidate for this could be 
the polar ferrimagnet GaFeC>3 . Recent studies have characterised the birefringence 
of this compound in the optical and x-ray ranges (Jung et al 2004; Kubota et al 
2004). Such studies suggest values larger than those induced by imposed fields in 
organometallic liquids: An ~ 10~ 4 (van Tiggelen et al 2005, van Tiggelen et al 
2006). This should amplify any vacuum effects if they are non-zero. We should 
point out that An rs A% drops with increasing temperature for GaFe03, vanishing 
with ferrimagnetic order at the transition temperature Tc ~ 225 K; further, at low 
temperatures small magnetic fields are required for a non-zero birefringence (Jung 
et al 2004). This means the radiometer should be placed in an evacuated chamber 
and maintained at temperatures and fields which maximise A\- Since light could 
cause the radiometer to turn, the chambers should be light-tight. Light from an LED 
laser placed on the top of each pane (see Fig. QJ>) would then allow a camera in the 
chamber to track any rotation of the radiometer. A mirror version of the radiometer 
should also be constructed and its rotation, or lack thereof, tracked. This allows 
to check the prediction that the rotation should reverse when the direction of the 
panes' optical axes is reversed, as well as allowing to detect unwanted drifts. 
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6. Discussion 

In this paper we argue for the experimental verification of Feigel's theory that 
the vacuum can transfer momentum to a fluid placed in strong crossed electric. 
The momentum transfer occurs because vacuum modes in such a fluid are no 
longer isotropic as in an ordinary dielectric. The crossed electric and magnetic 
fields Eextj B oxt change the symmetry of the dielectric fluid so that it behaves like 
a magnetoelectric medium with different refractive indices for waves propagating 
along or against the direction defined by E oxt x B oxt . Feigel's semi-classical argu- 
ment allows to derive the net momentum transfer caused by counterpropagating 
vacuum modes and the corresponding classical contribution due to relative motion 
in crossed fields. Our alternative derivation which confirms Feigel's result, once mi- 
nor inaccuracies in the original paper are corrected. Further, a new expression for 
the vacuum stress on the fluid is derived, predicting a Poiseuille flow in a tube, with 
maximum speed U max ~ 100/xm/s (2000 times larger than Feigel's original estimate 
of 50nm/s). This prediction contrasts with that of Obukhov & Hehl (2008) that a 
magnetoelectric slab in the vacuum experiences no net force. 

Two experiments are proposed to test the above predictions. In Feigel's original 
scenario of we predict fluid flow which can be measured from the flow rate from 
an organometallic fluid. Weaker flows could alternatively be measured from parti- 
cle tracking velocimetry with the same set-up. In the second experiment proposed, 
Obukhov and Hehl's null prediction for the force on a slab could be tested by mea- 
suring the rotation of a vacuum radiometer with panes made of a magnetoelectric 
material. 

Like all theories, Feigel's theory makes some assumptions which can and should 
questioned. The boldest assumptions, which we have adopted and are pivotal to 
our own realistic prediction, are the cut-off frequency assumption and the postu- 
late that vacuum modes see a magnetoelectric (assumptions 4 and 5, respectively, 
in Section |b|. The latter assumption amounts to assuming that the interaction of 
vacuum modes (virtual photons) with media is identical to that of light (photons). 
Since vacuum modes are, at least in part, electromagnetic, this seems reasonable; 
indeed, many QED effects both macroscopic (Casimir-Lifshitz forces) and micro- 
scopic (e.g. Lamb and Stark shifts) can be explained by considering the electro- 
magnetic interaction of matter with the vacuum (Milonni 1994). While microscopic 
approaches can be regularised using a cut-off based on the electron mass, as in the 
Lamb shift (Kawka et al 2010), Feigel's cut-off assumption has been criticised as 
improperly regularising the momentum integral, which, it claimed, makes the mo- 
mentum density incorrectly Lorentz variant. A proper dimensional regularisation, 
it has been argued, gives a null momentum density, go — (van Tiggelen et al 2005, 
van Tiggelen et al 2006). However, one could question the fact that the dimensional 
regularisation employed was pushed beyond the strict limits of its validity, yielding 
an erroneous result. When evaluating the surface tension of liquid helium films, 
Schwinger et al (1978) used a momentum cut-off based on interatomic distances 
and obtained predictions a factor of 3 larger than experiment. They stated a bet- 
ter agreement would have resulted by employing a microscopic model accounting 
for the short-distance physics not considered in their continuum model. We believe 
similar considerations apply to the Feigel effect. As to Lorentz invariance, it is the 
energy density spectrum, u (w), which needs to be Lorentz invariant (Boyer 1969). 
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Feigel theory preserves this invariance: since uq(u>) « go(w)c we see from (|2.6|) that 

A quantum microscopic model for the media would allow to properly account 
for absorption. As mentioned, high frequency modes will not sample a large enough 
number of atoms to experience a magnetoelectric response (which is the statistical 
product of many atoms). If this is the case, the vacuum contribution to momen- 
tum (pressure) from the absorption of high frequency modes is isotropic and can 
be neglected. If not, the exact consequences of absorption are hard to fathom in 
the absence of a microscopic model, but one would intuitively expect a smaller net 
transfer of momentum (from the inelastic collisions between the virtual photons 
and the medium, roughly speaking) . If the broadband spectrum of magnetoelectric 
(or effectively magnetoelectric) materials is similar to most dielectric/magnetic ma- 
terials, then the effect of absorption will be less important for smaller frequencies, 
which have smaller weight in the momentum density integral. By a similar argu- 
ment, the dispersion of the real part of the refractive index will matter more at high 
frequency. A model with a microscopically derived complex, frequency dependent 
refractive index will dispense of the need for a cut-off (or other regularisation). 

We should also comment on the agreement between Feigel's semi-classical deriva- 
tion and our own. The latter is also semi-classical in the sense that the form of the 
expectation value of the momentum density (our starting point) is implicitly given 
by quantisation of the classical momentum density. The question then arises as to 
what the appropriate momentum density should be. This is the so-called Abraham- 
Minkowski problem and, as many researchers have pointed out, it is not a problem at 
all, provided one is consistent about conserving total momentum (see, e.g., Loudon 
et al 2005). As can be seen from (|2.2p . Feigel uses the classical pseudo- momentum 
g = (DxB — Ex H)/47r. We directly quantise this, and evaluate (0|g|0) using: 

g = gAf - gA (6.1) 

where gM = (D x B + B x D)/47r is the Minkowski contribution, where photon 
momentum is proportional to refractive index n, and gA = (E x H + H x E)/47r, 
inversely proportional to it. Feigel on the other hand derives the difference in mo- 
mentum density for counterpropagating modes classically, and quantises the energy 
term which emerges. The two derivations thus differ only in the timing of the ap- 
plication of quantisation. Our (0|g|0) is taken to be proportional to n^x, which 
would suggest that we are using the Minkowski term only and not the expectation 
value of 16.11 Considering the expectation value of the Abraham form vanishes in a 
vacuum (0|g^|0) = (van Tiggelen et al 2005), we see that (0|g|0) = (0|g M |0), i.e. 
the momentum and its density scales like the refractive index, as we have assumed 
in our derivations. This Minkoswki pseudo-momentum has been shown to be the 
appropriate expression for the calculating momentum of matter caused by light 
(Peierls 1991). 

A comparison between the prediction by Obukhov & Hchl (2008) and that from 
Feigel's theory is also in order. The two scenarios differ in geometry. In Feigel's 
suggestion the fluid is unbounded, but in our revision this is true only in the di- 
rection along the optical axis. In Obukhov & Hehl's case, the magnetoelectric slab 
is bounded by its surfaces and it is the surface contributions which case the net 
vacuum stress to vanish. In our case, there are no surfaces so the vacuum stress is 
nonzero. In reality, there will be gradients on each side of the region of the fluid 
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where fields are applied, and it is conceivable that these also might cancel the vac- 
uum stress to vanish. The calculation of these contributions is beyond the scope 
of this paper, but would represent an interesting matter to pursue theoretically. In 
their paper Obukhov & Hchl state that the force on the slab by real photons (e.g. 
from counter propagating laser beams) should be nonzero. This seems to contradict 
their vacuum result and they do not make clear why real photons should behave 
different from virtual ones, though this may well be true. Both the experiments we 
suggest could be carried out with light and the expectation is that the fluid should 
move and the radiometer should turn when they are placed in counter propagating 
laser beams. The difference between real and virtual photons may be that the lat- 
ter are already involved in providing the magnetoelectric matter with its properties 
(fine structure etc.) so in some loose sense there are not enough virtual photons 
to cause a Feigel effect. Another objection might be that if the Feigel effect was 
real, we would be able to extract small amounts of energy from the vacuum with- 
out doing any work (e.g. using the radiometer or Feigel's wheels). It is not clear if 
second law of thermodynamics applies to the a system macroscopically coupled to 
the vacuum, since the latter is at zero temperature. On the other hand the third 
law states that absolute zero cannot be reached in a finite number of steps (Finn 
1993), so magnetoelectric media will always have a finite temperature. The vacuum 
momentum density and stress we have calculated may turn out to be be zero if 
finite temperature corrections are included. The Casimir effect does not however 
vanish at finite T (Lamoreaux 2005), so the Feigel effect could be could be simply 
another weird and wonderful consequence of quantum vacuum fluctuations. 

The Feigel effect has inspired many alternative theories, most adopting a regu- 
larisable Casimir geometry (van Tiggelen et al 2006, Birkeland & Brevik 2007). Any 
experimental predictions of these theories are beyond the reach of current instru- 
mentation with Feigel's original prediction recognised to be best candidate for an 
experimental test (van Tiggelen et al 2006). Given our realistic prediction is 2000 
times larger than the original, the case for an experimental test is even stronger. 
Testing Obukhov & Hehl's null prediction is also compelling. Both predictions in- 
volve a qualitative expectation which will hopefully be easy to detect: in one case 
the fluid moves (or not), in the other the radiometer stays still (or turns). Theo- 
ries can be challenged theoretically, but the final arbiter of theories in physics is 
experiment. This author is surprised that, some 7 years since it was proposed, an 
experiment to test Feigel's theory has yet to be carried out. We hope this work will 
stimulate experimentalists to find out if vacuum momentum can be asymmetrically 
transferred to matter. 
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Appendix A. Plane wave modes and refractive indices in a 
magnetoelectric with orthogonal symmetry 

In this appendix we derive expressions for the amplitudes and refractive indices of 
EM modes propagating in a magnetoelectric. The approach is very similar to the 
appendix of (Figotin & Vitebsky, 2001). When describing time-dependent fields in 
media, the following Maxwell equations are sufficient (see Chapter 41 of (Schwinger 
1998)): 



**.--!§ (AD 

(A2, 

In magnetoelectric materials with isotropic permittivity and permeability tensors 
(e = el; (x = fil) Equations (|A ip . (IA2[) are supplemented by the constitutive 
relations (Figotin & Vitebsky 2001, O'Dell 1970): 



D = eE + xH (A3) 

B = /iH + x T E (A 4) 

where the \ is the magnetoelectric susceptibility tensor, defined, in matrix form, 
by: 






Xxy 





Xyx 


















X = \ Xyx (A5) 
V / 

In (|A 15|) . Xxy, Xyx are magnetoelectric susceptibilities. They are a measure of the 
electrical polarisation caused by a magnetic field and the magnetisation caused by 
an electric fields applied along the x and y axes in the magnetoelectric. We assume 
plane wave solutions propagating in the z-direction: 

E kA (z, t) = ^oke i(fe ^- wt) e kA = £ kA e kA (A 6) 

with wave- vector k A , frequency w and polarisation vector e kA . For such waves (re- 
moving the suffices for clarity): 





' -Ey \ 




-l 





V x E = d z ( 


E x 


l-Mj 










J 












(A 7) 



Thus, defining: 




We can write: 



(A 8) 



VxE = (T<9 2 E (A 9) 
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Analogous considerations apply to V x H. Thus, since 9 2 E = ifcE and <9 4 B = — iujH 
and similarly for H and D, Equations (IA ip become: 



no-E = B (A 10) 

no-H = -D (All) 

where we have defined n = k{c/uj). To find the values of k for which the plane 
waves are solutions of Maxwell's equations. 

Substituting (|A~3|) and (|A~4l) into (IXlOl) and (fXTTj) and rearranging we find: 

H = — {na — x T )E (A 12) 

A* 

(na + x)H = -eE (A 13) 

Substituting for H from (|A 12[1 into (|A 13|) we find the following eigenvalue equation: 

AE = e/iE (A 14) 

where: 

(n — Xx ) 2 

X = -U,n+\)(na-\')= j " (n + Xxy ) 2 ) (Air,) 



(a) Mode Amplitudes and Indices of Refraction 

We look for linearly polarised plane wave solutions to Maxwell's equations. 
There are two independent polarisations. Consider first plane wave with E along 
ei = (1,0,0). The eigenvalue equation (IA 141) then entails: 

(ni - Xxyf \ / 1 \ / 1 \ 

(«i + Xxy? = en (A 16) 

/ \ / \ / 

So that the refractive indices for waves propagating in the k and — k directions with 
polarisation 1 are given by: 

«k : l = y/^+Xxy ( Al ?) 
"-k,l = -V^+Xxy ( Al 8) 

Analogously, for waves with E along e% = (0, 1,0) the eigenvalue equation (|A 141) 
requires: 

«k,2 = y/e+i-Xyx ( Al 9) 
«-k,2 = -^f^-Xyx (A 20) 

The magnetic field components corresponding to the above electric field modes are 
found using the Maxwell equation (|A 10p . For polarisation ei: 

/ o -m o \ / i \ / o 

Bi=nio-Ei= m } E! = m I E x (A 21) 
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and similarly for e 2 : 



B 2 = n2<j'Ei2 = 



-n 2 

J E 2 




(A 22) 



So finally, the four possible modes of propagation in a magnetoelectric are: 

[Eki,B k i] = -Eki[ei,n k ie 2 ], [E_ kl ,B_ k i] = -E-ki[ei,n_ k ie 2 ], (A 23) 
[E k2 ,B k2 ] = S k2 [e 2 , -n k2 ei], [E_ k2 ,B_ k2 ] = £- k2 [e 2 , -n_ k2 ei](A24) 

or, in Feigel's notation (l/E kX )(E x , E y , B x , B y ): 

(l,0,0,n kl ), (0,l,-n k2 ,0), 
(l,0,0,n_ k i), (0,l,-n_ k2 ,0). 



Appendix B. Lagrangian and Hamiltonian densities for 

magnetoelectric media 



In general the electromagnetic Lagrangian in ponderable media is given by: 

L = JJc-dr dt 
where the Lagrangian density is given by: 



C 



1 

4tt 



-E • D B • H 

2 2 



(Bl) 



(B2) 



A magnetoelectric material satisfies the constitutive relations (IA 3[) and (|A4p . These 
give D = D(E,H) and B = B(E,H). We can invert (fA~3l) to find H = H(E, B) 
and so D = D(E,B) from (TO)) : 



H = - (B x T E) 



D = eE 



-x(B-x T E 



(B3) 
(B4) 



Substituting (|B 3ft and (|B 41) into the Lagrangian density lB~2l we find: 



C = 



47T 



I E • ( e E + -xB - -xx T E J - ^B • (-B - - X T E 
2 \ [i (i J 2 V. ^ n 



or 



£ = — 

47T 



" 2^ + 2^ ' ^ + 2^ B ' ^ ~ 2^ E ' ^ TE > 



(B5) 



(B6) 



Letting Q = E ■ (£B) we see Q is a scalar (Q = Q T ) so that E • (£B) = B • (x T E) 
and, neglecting terms of order XxyiXj/ar 



£ = — 

47T 



»(X 2 ) 



Similarly, the Hamiltonian is given by: 



-ED+-BH 

2 2 



1 

4tt 



I eE 2 + — B 2 
2 2^ 



o( X 2 



(B7) 



(B« 
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quantity 


Feigel's value 


correct value 


magnetoelectric modes 


(1,0, 0,^) 
(1,0,0, -VP) 
(0,1, -v^t, 0) 
(0,1, v^S, 0) 


(l,0,0, v ^I + x« B ) 

(O.l.-Vep + Xw-.O) 
(0,1,VP + Xyx,0) 


momentum density of mode k 


2A X i(e+i)^ 




vacuum momentum density expectation 


32-/r :i "A m c 4 





Table 2. Trivially incorrect quantities in Feigel's derivation (Feigel 2004) an d the correct 

values derived in this paper. 

Appendix C. Inaccuracies in Feigel's original derivation 
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